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Abstract. We study the static behavior of a gas of atoms held in a one-dimensional 
lattice where two distinct electronically high-lying Rydberg states are simultaneously 
excited by laser light. We focus on a situation where interactions of van-der-Waals 
type take place only among atoms that are in the same Rydberg state. We analytically 
investigate at first the so-called classical limit of vanishing laser driving strength. We 
show that the system exhibits a surprisingly complex ground state structure with 
a sequence of compatible to incompatible transitions. The incompatibility between 
the species leads to mutual frustration, a feature which pertains also in the quantum 
regime. We perform an analytical and numerical investigation of these features and 
present an approximative description of the system in terms of a Rokhsar-Kivelson 
Hamiltonian which permits the analytical understanding of the frustration effects even 
beyond the classical limit. 
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1. Introduction 

Atoms in highly excited states — so-called Rydberg atoms — interact via power-law 
potentials, which in conjunction with an external laser drive give rise to intricate many- 
body phenomena. Recent experiments with Rydberg gases have revealed that the 
dynamical behavior of these systems is of surprising variety. Examples are the emergence 
of bistable behavior mg, the observation of correlated aggregation of excitations 
[3] and of coherent excitation transport [4, |5, [6!j. Strong interactions also become 
manifest in the structure of the ground state of Rydberg ensembles. Several theoretical 
works have predicted and studied the formation and the melting of crystalline phases 
d 0 0 EDI HD El El na ESI. But only recently it was shown that these crystalline 
states can be actually accessed experimentally [16]. 

Currently there is a surge of interest in atomic systems in which multiple Rydberg 
states are excited. One of the main motivations is the presence of an exchange interaction 
between Rydberg states of different atoms that results in coherent transport dynamics 
laiai, the non-local propagation of light na as well as a non-trivial collapse and 
revival dynamics [18, T9J. Interesting physics emerges also in the absence of exchange 
interactions. The main reason is that depending on the Rydberg states the interaction 
between two atoms can vary by one or more orders of magnitude [201121]. This feature 
is exploited for example in all-optical transistors for light pulses [22, [23] . 

In this work we explore the many-body ground states of an atomic lattice gas 
in which two Rydberg states — or species — are simultaneously excited. We focus 
on a situation where interactions are present only between atoms of the same species. 
This is reminiscent of the Potts model [24] , for which however only few studies exist that 
consider interactions that extend beyond nearest neighbors [25U26U2ZU2H! - We show that 
the ground state of the two-species Rydberg lattice gas features a surprisingly complex 
structure with a series of compatible to incompatible transitions. In the compatible 
case the two species can occupy the lattice such that they can both minimize their 
interaction energy independently. In the incompatible case this is not possible and it 
leads to mutual frustration. We provide analytical expressions for the regions of stability 
of the (in)compatible phases in the classical regime, i.e. in the limit of vanishing laser 
driving strength. The quantum regime is studied numerically as well as analytically with 
the help of an approximate Hamiltonian of Rokhsar-Kivelson form. These considerations 
show how frustration persists also in the quantum regime. 

2. Description of the system 

We consider a one-dimensional gas of atoms trapped in a lattice of spacing a with a 
single atom per site. The relevant internal level structure of the atoms is shown in Fig. 

The ground state | g) is coupled to two Rydberg nS-states (|s) and [?/;)) via two 
laser fields with Rabi frequencies and detuning A Q (a = s,w ). The fundamental 
interaction between two such atoms is to leading order given by the dipole-dipole 
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Figure 1. (Color online) (a) Relevant atomic levels. Two Rydberg nS-states |a), 
a = s,w (referred to as species) are excited from the ground state \g) by a laser with 
Rabi frequencies Q, a and detuning A a . Atoms in the state |s) (|'tc)) interact strongly 
(weakly). (b) Sketch of two minimum interaction-energy configurations for the strong 
and weak species at incompatible filling fractions, (c) Achieving a minimum energy 
disposition. The available lattice sites after placing the s-species form a distorted 
lattice with spacings of 1 and 2 times the original lattice spacing. Atoms of the w- 
species are placed in the available lattice sites according to the algorithm of Ref. 
[29] . (d) Examples of ground states with different filling fractions p a . Incompatible 
(compatible) cases are shaded in red (green). 


potential. However, for Rydberg atoms in nS-states, as considered here, this interaction 
results in a (second order) van-der-Waals energy shift [20J. For two atoms in the pair 
state \olol) = |a) |a) and separated by the distance R it reads V a (R) = Cq R~ 6 . The 
van-der-Waals coefficients Cjf 1 scale with the eleventh power of the principal quantum 
number n [30, [20l [31] . This generates strong interactions among Rydberg states and 
moreover permits to achieve a scenario in which the interaction between two atoms 
in the pair-state |ss) is much larger than the interaction between two atoms in \ww). 
Henceforth we focus on such a case, which is in practice achieved by choosing the 
principal quantum number of state |s) (the strong species) to be larger than that of state 
| w) (the weak species). Moreover, such large difference in the principle quantum numbers 
results in a strongly suppressed interspecies interactions [3.2] as experimentally shown 
in Ref. [33] . This is a rather generic feature of Rydberg atoms in the sense that one 
can find various combinations of suitable levels which satify such strong-weak hierarchy. 
Identification of these levels is atomic species dependent and the corresponding van-der- 
Waals coefficients can be calculated using perturbative techniques as described e.g. in 
[33] , To give an idea of the orders of magnitude of the Cq coefficients and how they 
depend on the principle quantum number, we provide a typical example in Fig. [2| where 
we compare the inter- and intra-species interactions for Rydberg states of rubidium. 

The Hamiltonian of the system can be written (within the rotating-wave 
approximation) as the sum of Hamiltonians of the strong and weak species, H = H w +H s , 
with 
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where n^ = |a) (a| and a x ^ = \a) (g\ + \g) (a\, and V a = Cq*> /a 6 . Note, that H w 
and H s do not commute due to the fact that both Rydberg species are excited from a 
common ground state. This is the origin of the compatible to incompatible transitions 
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Figure 2. (Color online) Comparison of the van-der-Waals CqCiq^u) coefficients 
describing the interaction strength between two rubidium atoms in Rydberg states with 
principal quantum number no and n respectively. Cq for the inter-species case is shown 
as a blue line (attractive interaction) and red circles (repulsive interaction). The black 
dashed line corresponds to the intra-species Cq, i.e. no = n. Green crosses represent 
a possible choice of Rydberg levels implementing the required hierarchy of interaction 
strengths assumed in our work. Identifying no = 65 with the strong species and n = 47 
with the weak one we get C^/Cq W ^ = 50.5 and /Cq SW ^ = 49, i.e. the hierarchy 
Cq S ^ Cq W ^ Cq SW ^ is well respected (we have denoted here by Cq SW ^ = Ce(65,47) 
the inter-species Cq coefficient). 


as well as the frustration effects which we will discuss in the following. 

3. Compatible to incompatible transition 

We start by studying the case — A a — 0 in the thermodynamic limit, i.e. an infinite 
lattice, where Eq. 0 reduces to the Hamiltonian of a classical two-species Potts-model 
with (convex) 1/r 6 interactions and inequivalent couplings ( V s ^ V w ). We are interested 
in finding the microscopic state which minimizes the interaction energy for given filling 
fractions of excitations p s and p w . For a single species and a general convex potential 
this problem was studied in [HE1 ESI ESI EZJ. Given the filling fraction the convexity 
of the potential forces the system into the most homogeneous configuration achievable 
accounting for the lattice constraint, as shown by Hubbard in Ref. [2Qj. 

When two species are present it is in general not possible to minimize at the same 
time the interaction energy of both. An example is given in Fig. [IJb) for p s = 1/3 
and p w = 1/2. Finding the actual ground state is simplified by our assumption 
Vs^Vyj. Here, the configuration of the s-species can be considered as “frozen” and not 
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constrained by that of the w-species[|J Atoms of the s-species will then arrange following 
the single species prescription of Ref. [29J, that is calling ri the distance between an 
s-atom and its Z-th nearest neighbor, the s-atoms will be distributed on the chain with 
distances satisfying r/ = [1/p s \ or \l/ p s ]. This set of constrains identifies a unique 
distribution for the atoms of the s-species, and will in general lead to a deformed lattice 
formed by the remaining empty sites on which the minimum-energy arrangement of the 
atoms of the w-species needs to be found. 

For the example shown in Fig. []Jb) the w-atom sitting in the doubly occupied site 
has to be moved to an empty site in a way that minimizes the increment in interaction 
energy. As shown in Fig. [IJc) the strong species leaves a distorted lattice with a density 
of empty sites given by 1 — p s = 2/3 and lattice spacings a and 2 a. The minimum 
interaction energy is then obtained by arranging the w;-atoms according to the single 
species prescription of Ref. [29] considering that the filling fraction of the w species on 
the distorted lattice is p w / (1 — p s ) = 3/4. This finally leads to the state depicted in 
Fig. B c). The mathematical proof of this method we provide in a separate publication, 
see Ref. [38]. Some examples of minimum interaction-energy configurations are shown 
in Fig. B d ) From these considerations one finds that two cases, which depend on the 
filling fractions p a , need to be distinguished: (i) The compatible case in which the two 
species assume their minimum-energy dispositions independently without interfering 
with each other, (ii) The incompatible case in which the strong species prevents the 
atoms of re-species from assuming their minimum energy configuration. 

Understanding whether two filling fractions are compatible is generally a hard 
task which requires the study of the full arrangement of excitations. In the following 
analytical analysis we will consider filling fractions of the form p a = 1 / q Q , where both q a 
are positive integers. Clearly, the minimum energy configuration is achieved when the 
a-atoms are arranged uniformly with distance q a . The filling fractions are compatible 
only if they share a common divisor (see [38j). In the incompatible case the distribution 
of the a;-atoms will be distorted and can be represented by repeated strings of repeated 
intervals of lengths ...q w (q w — l)(q w + 1) q w .... as shown for the example in Fig. [l](c) . 

Let us now turn to the discussion of the stability of the (in)compatible phases. So 
far we have assumed that the filling fractions p a of the individual species are externally 
imposed. However, in practice the filling fraction is controlled by the parameters A a in 
Hamiltonian ([Tj) which act as chemical potentials [HI Hi]. The question is then which 
state, or filling fraction, is actually stabilized in a given region of the A s — A w manifold. 
As the distribution of the s-species is effectively independent from that of the ic-atoms, 
we can study its stability with the single species method (HHUHTUH]. The stability of a 
given filling fraction p w on the other hand is less simple to analyze, as it depends on the 
configuration of the 5-atoms. We can obtain an analytical result for the transition to an 
incompatible state in the thermodynamic limit starting from a compatible arrangement 

f This is strictly speaking true only in the asymptotic limit V s /V w oo. Otherwise an additional 
assumption has to be made on the filling fractions, namely V s /V w 1 and p s /p w > 1? the situation 
we consider here. See also [38] for details. 
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Figure 3. (Color online) Exact diagonalization results, (a) Density plot showing the 
transverse magnetization of the strong species m s in a parameter regime A s where the 
filling fraction in the classical limit p^ 1 , i.e. at Q s = 0, is 1/2. The white contours 
show lines of equal p s . (b) Frustration energy ©> plotted along A s = 102, i.e. the 
p s = 1/2 contour of (a). The blue and red lines correspond respectively to a compatible 
(p^ 1 = 1/4) and an incompatible (p^ = 3/8) classical filling fraction. The value of p w 
is indicated explicitly in the figure for some plateaus. The inset magnifies the first 
small jump in p w in the incompatible cases. (c,d) Filling fractions of the weak species 
for compatible [p^ 1 = 1/4, (c)] and incompatible [p^ 1 =3/8, (d)] case. The data were 
obtained with V s = 10014, = 100. 


of atoms for which v = q w /q s is integer (see Section 6 in [38]). The region of stability is 
delimited by 

A^ 1 =F l(q w — v) l(q w — 2u) lis 

V w ~ ^ [ ( lq w ) 6 (lq w T l) 6 (lq w T 2) 6 

On the one hand a transition to an incompatible state takes place when A w = , the 

value of A^ for which the introduction of one more excitation is energetically favorable 
when keeping A s fixed. In this case the distribution of the Z-th nearest neighbor rc-atoms 
which in the compatible phase is homogeneous, e.g. atoms at distance lq w , is deformed 
by the introduction of Iv distances q w — 2, and l(q w — 2v) distances q w — 1. On the other 
hand there is a transition at A w — A^ at which the w species loses an excitation and 
the compatible distribution is deformed by the introduction of Iv distances q w + 2 , and 
l(q w — 2u) distances q w + 1. 


( 2 ) 
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4. Quantum fluctuations and frustration energy. 

We now address the question as to how quantum fluctuations introduced by a laser of 
finite driving strength affect the compatible to incompatible transitions and in particular 
the emergence of frustrated states. To this end we consider a non-zero value of Ll s . This 
changes the state of the s-species from a classical one to a superposition of configurations 
with different excitation number. We are interested in how this impacts the classical 
arrangements of the re-species, i.e. at Q w = 0. To quantify this we introduce the 
frustration energy 

E,r = 47(A.,*; Si,) - 4p (AJ, (3) 

which measures by how much the strong species prevents the w-species from reaching 
its minimum energy configuration if it were alone. Here E^sp (^s,w', fl s ) = (H w ), i.e. 
the expectation value of H w [see Eq. |lj)] and (A u ,) is the energy of the classical 
configuration of the re-species in the absence of s-atoms. 

In the classical limit, fl a = 0, the frustration energy is zero when the detunings 
A q are chosen such that they stabilize filling fractions p a which are compatible. iq r 
becomes in general larger with increasing fl s as increasing density fluctuations in the 
s-species force the re-atoms to assume configurations with increased energy. In order to 
study this behaviour in more detail we diagonalize Hamiltonian (JT|) for a chain of length 
L = 8 with periodic boundary conditions. This approach has the typical drawbacks of 
a small scale numerical study: the presence of finite size effects (though minimized by 
the periodic boundaries), and the fact that one is limited to filling fractions of the form 
p/L for p < L. We will see, however, that the results provide an intuitive understanding 
of the physics at work. Note furthermore, that the considered relatively small system 
size in fact comes close to what is currently realizable experimentally in the context of 
Rydberg atoms mm- 

For our numerical study we focus on a regime where the strong species in the 
classical limit forms a crystal with filling fraction pf = 1/2. With increasing Q s this 
crystal melts as is seen in Fig. |3](a) where we show a density plot of the transverse 
magnetization m s = (1/L) (oy' s) ). The magnetization displays the typical lobe- 

structure H, and the formation of a (longitudinally) paramagnetic state (m s = — 1) 
at large fl s . Here it is evident that the state of the strong species is formed by a 
superposition of states with different number of excitations. 

Let us now have a look at the frustration energy when the weak species is added. 
The corresponding data is shown in Fig. |3](b). In the classical limit (Q s = 0) we choose 
A s = 102 which corresponds to half filling of the strong species, pf = 1/2, and consider 
values of A u , which stabilize the compatible case pf = 1/4 (A u , ~ 0.0159 blue line) 
and the incompatible case pf = 3/8 (A„, ~ 0.0314, red line), respectively. We first 
determine the values of A u , numerically as the midpoints of the stability regions. In 
the compatible case one can use Eq. 0 to determine A w . In the incompatible case 
on the other hand Eq. (|2j) cannot be applied directly. In the finite size case though, 
as previously mentioned, only fillings of the form p/L are accessible, and as such the 
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region of stability of p^ = 3/8 extends from the upper boundary A^ for p w = 1/4 to 
the lower boundary for pff = 1/2. These cases are compatible such that we can 
use again Eq. (|2j) to estimate the midpoint of this region of stability. We find that 
the values obtained analytically from Eq. 0 are in remarkable agreement with the 
numerically determined values despite the small size of the system. At = 0, lq r is 
zero for compatible densities and finite for incompatible ones. With increasing f2 s , the 
w-species becomes more and more frustrated which is reflected in an increase of E R . 
Interestingly, this increase is step-wise and each step is accompanied by a change of the 
filling fraction p w of w-atoms, whose value is provided in Fig. [dj^b) for each plateau of 
E{ t . 

Furthermore, we show the behavior of the weak filling fraction in the entire fl s — A s 
plane in Fig. [3|c,d). Note, that contrary to p s [shown as contours in Fig. la)], ft 
exhibits a lobe like structure which is not symmetric. This can be understood as follows: 
For a given finite the filling fraction p s decreases with decreasing A s . The resulting 
smaller number of s -atoms permits the accommodation of a larger number of w-at orris. 
This effect is stronger the larger Q s leading to an increasing asymmetry of the lobes 
in Figs. [3](c,d) . Moreover, for any A s , p w decreases stepwise with increasing fl s and 
eventually vanishes. The decrease in p w can be qualitatively explained as follows: As 
Q s increases the ground state begins to contain basis states with s-atom numbers that 
are larger than in the classical case. This forces the w-atoms to assume a lower filling 
fraction. The fraction of admixed basis states with increased s-atom number becomes 
larger the larger A s [see contours in Fig. [3j[a)] which results in the asymmetry of the 
lobes in Figs. 0C, d). 

We expect the qualitative features, namely the increase of frustration with f l 3 to 
hold also in the thermodynamic limit. To study this a semi-analytical treatment of 
quantum fluctuations based on a perturbative expansion of the ground state around the 
classical limit could in principle be performed. This task is, however, quite involved 
due to the fact that the filling fraction of the weak species is in fact a function of Q s . 
Nevertheless, one can still gain analytical insight into the quantum regime by considering 
the description of the system in terms of an approximate Hamiltonian, as is shown in 
the following. 

5. Rokhsar Kivelson approximation 

In the following we conduct an approximate analytical study of the ground state by 
means of a Rokhsar-Kivelson (RK) Hamiltonian Hrk- This generalizes the approach 
used in Refs. mmm for characterizing the statics of a Rydberg gas. The central 
idea is to approximate the Hamiltonian 0 as a sum of local positive semi-definite 
projective Hamiltonians H RK = Ylk=i ^k- In this case a state |G) can be found which is 
annihilated by each of the local Hamiltonians, i.e. hk\G) = 0, and thus represents the 
ground state of H RK . Under the RK approximation the ground state |G) is given by 
a superposition of classical configurations of a gas of hard-core polymers. As is shown 
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Figure 4. (Color online) (a) RK manifold of the s-atoms (solid blue line) and the 
re-atoms (red dashed line) for V s = 10014; = 100. The RK manifold is parameterized 
by the parameters £ a which increase when going from right to left along the manifolds. 
The inset show a magnification of the main figure in the region of small Q a . (b) Filling 
fractions p a = (1 /L)Ywk=i (G| \ G) as function of £ w in the absence of s-atoms. 

Here p w saturates to its unfrustrated value 1/2. The blue circles and red triangles 
represent the filling fractions obtained by diagonalizing exactly Hamiltonian 0 for 
L=6 with periodic boundary conditions. We chose V s = 10014, = 100, and used the 
definitions of £ s?w , to fix and Eq. (|6| to fix A SjW . (c) In the presence of s- 

atoms the re-atoms can only assume a maximum filling fraction p s = 1/3 and form a 
frustrated state with exponentially many configurations. The inset shows two example 
configurations with filling p w = 1/3. Here red ellipsoids (blue rectangles) sketch the 
dimers (trimers) of the RK construction. The blue circles and red triangles are the 
exact results drawn for comparison, and are obtained with the same method explained 
for panel (b). 


below, this procedure neglects the long-range tails in the interactions. As a consequence 
the densities that the different species assume in the classical limit need to be adjusted 
by hand. 

We consider a situation where in the classical limit the densities of the weak and 
strong species are given by p w = 1/2 and p s = 1/3. In this case the we can think of 
atoms in terms of polymers extending over sites at maximal distance R w — 1 (dimers) 
and R s — 2 (trimers), respectively (see Ref. |4TJ for detail). The polymer analogy 
is equivalent to the hard-core conditions — 0 and n^n^!± 1 = n^n j^ 2 = 0. 

Neglecting interactions on distances larger than R a + 1 and making the hard-core 
conditions explicit also in the laser-excitation part of the Hamiltonian HUSH , we can 
rewrite Hamiltonian © as 

L 

h= Y, {- Q ^a + a4 Q) + [ A « - - a*(2 R a + qej d Q) 

k=l 1 a=s 1 w 


(a) (a) 

n k n k+R +1 
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where the local Hamiltonians are given by 


4“’ = n t 


(«) 


c(a) 




(<?) 


+ L ln i a) 


(5) 


The //|, 0 ) are the above-mentioned positive semi-definite operators which are proportional 
to projection operators and whose sum yields the RK Hamiltonian. The projection 
operators implement the hard-core condition by ensuring that dimers (trimers) 
do not occupy neighboring (next-neighboring) sites: 11^ v) = h-ihti and n£ } = 

(s) 0 ) 0 ) 0 ) (a) . (a) 

Pk-2Pk-lPk+lPk+2 Wlth Pk = 1 - U k • 

Hamiltonian Q is not yet of RK form since there are a number of additional terms. 
For specific parameter choices, however, these terms vanish up to trivial constants. 
In order to make this more transparent we have introduced (following Ref. [TOj) the 
parameters f a . Those can be chosen freely since the Hamiltonian does not depend on 
them, as can be seen when multiplying out all the individual terms. 

The RK form is assumed when choosing £ Q = f a V a /Q a with f s = 1/(2 x 3 6 ) and 
f w = 1/2 6 and requiring the parameters f2 a , A a , V a to accomplish 

2 A 


a 


+ 9a > 


( 6 ) 


JaV a J f a V c 

where g s = 5 and g w = 3. The first condition eliminates the two terms in the second 
line of Eq. 0 while the second one leads to a cancellation of the terms proportional to 
A a and il a . Both conditions define the so-called RK-manifold [see Fig. [4](a)] on which 
Hamiltonian (|4]) assumes the RK form 


Hrk (f^s^s T L T ^ ' £l a h 


(a) 


(7) 


k=l : cx=s,w 

The ground state | G) of H R k is written as 


N 


k =1 


\G) =-[[(!-WPcC-W 


ds)(s)t 
L k c k 


| 0 ), 


( 8 ) 


where |0) = <S>iLi \9)i an d = |a) (g |. The parameters can be then thought of as 
the fugacities of the dimers and trimers gases, such that large (small) values result in 
a ground state with high (low) density of the respective species. The ground state 
can be represented explicitly as matrix product state (see e.g. HU)- We can now use 
the ground state Hrk to get an idea of the effect of quantum fluctuations on crystalline 
states. We find that frustrated states emerge also away from the classical limit, i.e. for 
non-zero Cl a . 

By construction Hrk is frustration free, i.e. the frustration energy © is always zero. 
In fact this is an artifact stemming from the omission of long-range tails. Nevertheless, 
the interspecies frustration still manifests in the behaviour of the filling fractions, Fig. [4j 
In Fig. |4](b), where we set < 0 in order to have a low density of dimers, we find that 
indeed only dimers are occupying the lattice and that they fill it completely in the limit 
—> oo. This corresponds to a unique (translationally invariant) crystalline ground 
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state with p w = 1/2. On the other hand in Fig. |4]^c) we fix the parameter £ s such that 
the strong species (trimers) fills up the lattice at density p s = 1/3. We now increase 
£,,, which increases the density of the weak species. However, instead of a saturation 
at density 1/2 — as in the absence of the strong species — we find that p w tends to a 
value of 1/3. This is a signature of inter-species frustration. The weak species does not 
crystallize but is in a superposition of exponentially many states with density 1/3, two 
of which are shown in the inset of Fig. [4|h). 

6. Summary and outlook 

We have investigated the statics of a two component Rydberg gas at zero temperature, 
with particular focus on the emergence of frustration effects. In the classical regime we 
identified regions in parameter space where the two Rydberg species form compatible 
and incompatible arrangements. In the quantum case we introduced the frustration 
energy to quantify the degree of frustration. We found that this quantity shows a 
staircase pattern as the system visits different phases characterized by different filling 
fractions. Finally, we performed an approximate analytical study of the quantum regime 
by means of an RK Hamiltonian, which further corroborated the existence of frustrated 
states. 

Our analysis provides insights into the complexity of ground state structures in 
multi-component Rydberg gases. In the future it would be interesting to extend the 
work to situations without the strong-weak hierarchy that was assumed here. Moreover, 
an extension of the analysis to Rydberg gases with more than two relevant excited states 
and/or to high-dimensional lattices would be desirable due to the number of Rydberg 
lattice experiments that recently has become available. 
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